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LONG-TERM  GOALS 

Achieving  reliable  communications  with  high  data  rates  over  underwater  acoustic  communication 
(UAC)  channels  has  long  been  recognized  as  a  challenging  problem  owing  to  the  scarce  bandwidth 
available  and  the  time-varying  channels  with  large  delay  spread.  Fortunately,  the  underwater 
environment  with  multipath  scattering  makes  it  possible  for  multi-input  multi-output  (MIMO)  spatial 
multiplexing  schemes  to  achieve  high  data  rate  UAC.  The  major  challenges  of  the  high  data  rate 
MIMO  UAC  scheme  include  time-varying  channels,  long  delay  spreads,  mutual  interferences  among 
the  data  streams  transmitted  by  the  different  transducers,  and  Doppler  scaling  effects  due  to  relative 
movement  between  transmitter  and  receiver  platforms.  These  challenges  require  short  and  effective 
training  sequences,  reliable  temporal  resampling  schemes,  efficient  channel  estimation  and  tracking 
algorithms,  and  innovative  equalization  and  decoding  strategies.  Our  long-term  goals  are  to  efficiently 
and  effectively  tackle  these  challenges. 


OBJECTIVES 

The  objectives  of  our  proposed  research  are:  1)  to  synthesize  effective  and  efficient  training  sequences 
with  cyclic  prefixes  for  the  accurate  estimation  and  prediction  of  UAC  channels,  2)  to  develop  a 
temporal  resampling  scheme  that  can  sufficiently  compensate  for  the  Doppler  scaling  effects,  3)  to 
devise  robust  and  computationally  efficient  sparse  channel  estimation  and  tracking  algorithms  that  can 
be  used  to  provide  statistical  properties  of  the  estimated  UAC  channels,  4)  to  design  efficient  detection 
scheme  that  can  effectively  perform  interference  cancellation  in  a  MIMO  system,  and  to  provide 
channel  updating,  equalization,  and  decoding  sehemes  with  good  initial  conditions  to  allow  improved 
turbo  processing,  and  5)  to  continue  the  ongoing  in-water  experimentation  and  algorithm  evaluation 
activities  by  continuing  our  collaborations  with  the  Woods  Hole  Oceanographic  Institution  (WHOI) 
and  others. 

APPROACH 

Training  sequence  design:  We  have  recently  introduced  a  new  cyclic  algorithm  (called  CAN  =  CA- 
Ncw)  for  the  local  minimization  of  the  integrated  sidelobe  level  (I SL)  metric  for  efficient  aperiodic 
sequence  designs.  CAN  is  based  on  fast  Fourier  transform  (FFT)  operations  and  can  be  used  virtually 
for  any  practically  relevant  values  of  the  sequence  length,  which  can  be  up  to  a  million  or  even  larger. 
We  have  also  extended  CAN  to  generate  perfect  periodic  sequences,  which  are  the  so-called  constant 
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amplitude  zero  auto-correlation  (CAZAC)  sequences.  This  CAZAC  sequence  design  problem  goes 
back  to  Norbert  Wiener  or  even  Gauss!  The  mathematicians  have  shown  that  for  some  sequence 
lengths,  there  are  an  infinite  number  of  perfect  periodic  sequences.  We  have  solved  a  long-standing 
perfect  periodic  sequence  design  problem  -  finding  such  sequences  (including  long  sequences) 
efficiently  and  effectively.  We  solved  it  via  exploiting  FFTs. 

The  extended  CAN  algorithm  is  referred  to  as  periodic  CAN  (PeCAN).  Unlike  most  existing  sequence 
construction  methods  which  are  algebraic  and  deterministic  in  nature,  we  start  the  iteration  of  PeCAN 
from  random  phase  initializations  and  then  proceed  to  cyclically  minimize  the  desired  metric.  In  this 
way,  through  different  random  phase  initializations,  we  can  find  as  many  different  waveforms  as  we 
may  ever  want!  The  so-generated  sequences  are  difficult  to  guess  by  the  foe  and  hence  are  especially 
useful  as  training  sequences  or  as  spreading  sequences  for  covert  UAC  applications.  We  will  use 
PeCAN  sequences  for  more  in-water  experimentations  to  demonstrate  their  effectiveness. 

Temporal  Resampling:  In  mobile  MIMO  UAC,  the  presence  of  Doppler  effects,  owing  to  the  relative 
motions  between  the  transmitter  and  receiver  platforms,  induces  temporal  scaling  (stretching  or 
compression)  to  the  transmitted  signals.  The  Doppler  scaling  effects  can  be  sufficiently  compensated 
for  by  converting  the  original  acoustic  channel  subject  to  both  inter-symbol  interference  (IS1)  and 
Doppler  scaling  effects  to  an  ISI  channel  via  resampling  the  received  measurements.  Although  the 
Doppler  scaling  effects  can  be  largely  mitigated  after  such  a  temporal  resampling  process,  the  residual 
Doppler  still  causes  frequency  shift  on  the  received  measurements,  which  requires  the  employment  of 
channel  estimation  algorithm  that  can  accurately  estimate  the  resulting  frequency  modulated  acoustic 
channel. 

Channel  Estimation:  Through  simultaneous  transmissions  using  multiple  transmitters,  MIMO 
systems  offer  enhanced  reliability  and/or  increased  data  rates  compared  to  their  single-input 
counterparts.  However,  in  MIMO  systems,  the  multiple  transmissions  lead  to  mutual  interferences  at 
the  receiver.  To  successfully  recover  the  transmitted  data  sequences,  accurate  channel  estimation 
techniques  play  crucial  roles.  We  focus  on  MIMO  UAC  over  sparse  acoustic  channels  suffering  from 
1S1  and  frequency  modulations.  We  first  present  an  extension  of  our  recently  introduced  SLIM 
algorithm,  which  stands  for  Sparse  Learning  via  Iterative  Minimization,  to  estimate  the  sparse  and 
frequency  modulated  acoustic  channels.  The  extended  algorithm  is  referred  to  as  generalization  of 
SLIM  (GoSLIM).  GoSLIM  assumes  that  at  each  receiver,  the  channel  taps  for  all  the  transmitters 
experience  the  same  Doppler  frequency,  but  different  receivers  experience  different  Doppler  shifts.  We 
further  simplify  this  channel  model  by  assuming  that  the  channel  taps  for  all  the  transmitter  and 
receiver  pairs  experience  the  same  Doppler  frequency.  Accordingly,  a  variation  of  GoSLIM,  referred 
to  as  GoSLIM-V  (V  stands  for  variation),  is  presented.  We  show  that  compared  to  GoSLIM,  GoSLIM- 
V  not  only  reduces  the  overall  complexity  involved  in  the  channel  estimation  stage,  but  also  slightly 
improves  the  detection  performance.  This  comparison  between  GoSLIM  and  GoSLIM-V  suggests  that 
proper  channel  modeling  is  critical  for  effective  and  efficient  mobile  MIMO  UAC. 

Receiver  Design:  The  channel-induced  phase  shift  should  be  first  compensated  for  using  the  Doppler 
frequency  estimate.  Such  phase  compensation  task,  along  with  the  aforementioned  temporal 
resampling  process,  effectively  converts  an  acoustic  channel  subject  to  both  ISI  and  Doppler  scaling 
effects  to  an  ISI  channel,  which  allows  for  the  employment  of  various  equalization  techniques  that  can 
effectively  combat  ISI.  We  use  a  linear  minimum  mean-squared  error  (LMMSE)  based  filter  for 
symbol  detection.  In  a  MIMO  setup,  on  top  of  ISI,  multiple  simultaneously  transmitted  signals  act  as 
interferences  to  one  another.  Therefore,  interference  cancellation  scheme  also  plays  a  critical  role  in 


the  overall  detection  performance.  A  hard  decision  based  interference  cancellation  scheme,  including 
vertical  BLAST  (V-BLAST)  and  RELAX-BLAST,  subtracts  the  hard  decisions  of  detected  signals  out 
from  the  received  measurements  to  aid  the  detection  of  the  remaining  signals.  By  combining  V- 
BLAST  with  the  cyclic  principle  of  the  RELAX  algorithm,  RELAX-BLAST  provides  superior 
detection  performance  over  V-BLAST  at  the  cost  of  slightly  increased  complexities. 


The  detection  performance  can  be  further  enhanced  by  employing  a  soft  interference  cancellation 
scheme,  including  Turbo  equalization.  For  a  receiver  employing  Turbo  equalization,  both  the  equalizer 
and  decoder  involved  are  configured  as  soft-input  soft-output.  The  detection  performance  improves  as 
the  soft  information  cycles  between  the  equalizer  and  decoder.  The  main  drawback  of  the  Turbo 
equalization  scheme  is  the  increased  computational  complexity  compared  to  its  hard  decision  based 
counterparts.  To  address  this  problem,  we  consider  a  low  complexity  approximation  of  soft-input  soft- 
output  equalizer,  which  enjoys  a  computational  complexity  comparable  to  RELAX-BLAST  and 
provides  only  slightly  degraded  detection  performance  compared  to  a  directly  implemented  equalizer. 


WORK  COMPLETED 


The  received  measurements  should  be  first  temporally  resampled  to  mitigate  the  Doppler  scaling 
effects.  We  take  advantage  of  the  preamble  and  the  postamble  of  a  data  packet  to  estimate  the  Doppler 
scaling  factor  via  channel  estimation.  Specifically,  by  using  a  channel  estimation  algorithm  for 
acoustic  channels  suffering  from  Doppler  scaling  effects,  we  can  obtain  the  channel  impulse  responses 
(CIRs)  from  the  two  measurement  segments  in  response  to  the  preamble  and  postamble.  The  time 
duration  change  Td  imposed  on  the  packet  can  be  inferred  from  the  tap  shift  of  the  principal  arrivals  of 


»  T  +T. 

these  two  CIRs.  Then  the  Doppler  scaling  factor  estimate  is  computed  as  J3  =  — - — ,  where  Ta 


represents  the  duration  of  the  same  packet  at  the  transmitter  side.  Since  (i  can  never  be  perfectly 
accurate,  after  temporal  resampling,  Doppler  frequency  shifts  still  exist.  Accordingly,  the  resulting 
acoustic  channels  can  be  reasonably  modeled  as  frequency  modulated  channels. 


By  extending  the  SLIM  algorithm  to  deal  with  Doppler  frequency  modulation,  we  have  devised  the 
GoSLIM  algorithm  to  estimate  the  underlying  CIRs  and  Doppler  frequency  in  a  joint  manner.  Before 
proceeding  to  discuss  GoSLIM,  it  is  instructive  to  review  a  simpler  case,  i.e.,  the  channel  estimation 
problem  in  stationary  channels  using  the  SLIM  algorithm.  We  consider  a  M1MO  UAC  system 

equipped  with  N  transmitters  and  M  receivers.  For  stationary  frequency-selective  channels,  the 

v 

channel  estimation  problem  at  the  wth  receiver  is  given  by  ym  =  V X„hnm  +em ,  where  ym  is  the 

n= 1 

received  measurement  vector,  Xn  is  constructed  using  the  training  sequences  in  the  training-directed 
mode  or  the  previously  detected  symbols  in  the  decision-directed  mode,  hnm=[hnmr...,hnmli]T 
characterizes  the  C1R  of  length  of  R  between  the  h01  transmitter  and  the  m'h  receiver,  and  em  is  the 
additive  noise.  The  channel  estimation  problem  can  be  expressed  more  compactly  as  ym  =Xhm  +em , 
where  X  =  [X,,...,X v] and  hm  =  [h'm,...,h'Vm]'  .  We  consider  SLIM  for  estimating  hm  given  ym  and 
X  and  assume  that  the  noise  vector  em  contains  a  circularly  symmetric  independent  and  identically 


distributed  (i.i.d.)  complex-valued  Gaussian  random  variables  with  zero-mean  and  variance  rjm .  Let 
Pn.m.M  the  variance  of  hnmr,  and  define  P and  pm  =  [p;,m,-.p7v,Jr  •  Define 
Pm  =diag(pm)  as  the  covariance  matrix  of  the  channel  taps  seen  by  the  in"1  receiver  and  \nmr  as  the 
column  of  Xthat  corresponds  to  hnmr.  SLIM  is  user  parameter  free,  making  it  easy  to  use  in  practice. 
The  SLIM  algorithm  can  be  implemented  as: 

x"  v 

•  Step  0:  Initialize hnm=—r—:L —  . 

r  n,m,r  // 

XMm  X  _ 
ntm,r  ntm,r 


Step  1 :  With  the  most  recent  h„ ,  calculate  p„  „  =  I h„  m  rf . 

i  m  ’  /  n,m,r  \  n,m,r  | 

Step  2:  With  the  most  recent  Pm  and  //m,  calculate  hm  =[XllX  +  ^mPm~l )  X"ym  . 

IL-xhil2 

Step  3:  With  the  most  recent  hm,  update,  ijm  =  - - — ,  where  d  denotes  the  length  of  ym. 

d , 


•  Iterate  Steps  1-3  until  convergence. 


For  frequency-selective  channels  further  modulated  by  a  Doppler  frequency  fm  at  the  m'h  receiver,  the 
channel  estimation  problem  becomes  ym=AmXhm+em.  The  vectors  ym ,  hm ,  and  em  have  already 
been  defined  previously  and  Am  =  diag [\,e~2j*fmT‘ ,•  •  •  represents  the  Doppler  shift  matrix. 

In  Am ,  fm  is  the  Doppler  frequency  to  be  estimated  and  Ts  represents  the  symbol  period.  Given  ym 

and  X ,  we  consider  the  extension  of  the  SLIM  algorithm,  referred  to  as  generalization  of  SLIM 
(GoSLIM)  to  estimate  the  CIR  hm  and  the  underlying  Doppler  frequency  fm  jointly.  The  GoSLIM 
algorithm  can  be  implemented  as: 

•  Step  0:  Initialize  h  =  and /„  =  0 . 

X  X 

n,m,r  n,m,r 


9 


Step  1 :  With  the  most  recent  hm ,  calculate pnmr  =\hn„  r|2 . 

Step  2:  With  the  most  recent  Am ,  Pm ,  and  tjm ,  calculate  hm  =  (x"X  +  7/mPm_l )  ‘  ( AmX)"  ym . 

Step  3:  With  the  most  recent  hm ,  we  calculate  xm  =  Xhm .  Denote  xm{i )  and  ym(i)  as  the  ilh 
element  of  xm  and  yra  ,  respectively.  Then,  the  Doppler  frequency  can  be  calculated  as 

(  Jy  \ 


fm  =  arg  max  Re 


-2jnfmT,(i-\) 


,  where  zm(i)  =  y‘m(i)xji) . 


/, 


J 


Step  4:  With  the  most  recent  Am  and  hm ,  calculate  tjm  = 


|ym-AmXhjr 


•  Iterate  Steps  1-4  until  convergence. 


The  channel  estimation  problem  ym  =AmXhm  +em  corresponds  to  an  assumption  that  the  NR  channel 
taps  seen  by  the  mlh  receiver  experience  the  same  Doppler  frequency  fm ,  but  the  frequency  value 
could  vary  at  different  receive  hydrophones.  We  consider  a  further  simplified  assumption  that  the 


Doppler  frequency  is  the  same  across  M  receivers,  i.e.,  /  =  /,  =■■■■  =  fM.  Accordingly,  replacing  fm 
in  the  Doppler  shift  matrix  Am  by  /  yields  a  matrix  A  =  diag^l,e_2y'r/r,,-'-,e’:"ry7’(‘,,~l)j  that  is 


independent  of  the  receiver  index  m  .  Accordingly,  the  measurements  at  the  m"'  receiver  can  be 
expressed  as  ym  =  AXhm  +  em .  Stacking  the  measurements  from  all  the  receivers  follows  y  =  AXh  +  e , 

where  y  =[y[,...,y[/]'  ,  h  =  [hf,...,hj,]r  ,  e  =  [ej‘,..,ej/f  ,  A  =  lMxM®A,  and  X  =  1UM  ®  X  (I 

denotes  an  identiy  matrix  and  ®  denotes  the  Kronecker  product).  Again,  we  assume  that  e  contains 
circularly  symmetric  i.i.d.  complex-valued  Gaussian  random  variables  with  zero-mean  and  variance  . 
Then,  the  variation  of  GoSLIM,  referred  to  as  GoSLIM-V.  aims  to  estimate  /  and  h  given  y  and  X  . 
Note  that  unlike  GoSLIM.  which  can  be  employed  at  each  receiver  to  estimate  the  channel  in  parallel. 
GoSLIM-V  implicitly  suggests  that  the  measurements  acquired  at  different  receivers  should  be 
assembled  in  a  central  processor  before  performing  channel  estimation.  Moreover,  GoSLM-V 
simultaneously  estimates  the  CIRs  among  all  of  the  MN  transmitter  and  receiver  pairs.  The  GoSLIM- 
V  algorithm  can  be  implemented  as: 

x11  v 

•  Step  0:  Initialize  hn m  r  =  —  and  /  =  0  . 

’  ’  v  X 

n,m,r  ntmtr 


Step  1 :  With  the  most  recent  h ,  calculate p„mr=  | hn m r 
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Step  2:  With  the  most  recent  A ,  P  and  ^.calculate  h=^XHX  +  ^P  (AX)"  y  . 


Step  3:  With  the  most  recent  we  calculate  x„  =  Xhm.  Denote  xm(i)  and  y„(i)  as  the  ilh 

element  of  xm  and  y,„  ,  respectively.  Then,  the  Doppler  frequency  can  be  calculated  as 

A  /  m  \ 


/  =  arg  max  Re 

/ 


U  EZm(0 


/=!  V  m=\ 


where  zm{i)  =  yji)xm(i) . 


Step  4:  With  the  most  recent  A  and  h,  calculate  rj  = 


jy-AXh 

Mdy 


Iterate  Steps  1-4  until  convergence. 


Empirical  experience  indicates  that  the  computational  bottleneck  of  both  GoSLIM  and  GoSLIM-V  is 
in  the  update  of  the  Doppler  frequency.  Since  GoSLIM  is  employed  at  each  receiver  to  conduct 
Doppler  frequency  estimation  in  parallel  while  GoSLIM-V  estimates  the  common  Doppler  frequency 
across  M  receivers,  M  times  complexity  reduction  is  achieved  by  employing  GoSLIM-V  when 
updating  the  Doppler  frequency  estimation.  As  a  consequence,  GoSLIM-V  is  computationally  much 
more  efficient  than  GoSLIM  in  a  M1MO  configuration,  especially  with  a  large  number  of  receive 
elements.  Due  to  the  reduced  number  of  unknowns,  the  GoSLIM-V  data  model  is  more  parsimonious 
than  that  of  GoSLIM,  which  could  enhance  the  symbol  detection  performance  if  the  model  is 
reasonably  accurate. 


For  Doppler  modulated  frequency-selective  channels,  once  the  estimate  of  the  Doppler  frequencies 
{/m)m= i  (or  Doppler  frequency  /)  is  available,  the  receiver  can  compensate  out  the  adverse  effects 

induced  by  the  Doppler  shift  and  effectively  convert  the  original  frequency-modulated  channels  to 
stationary  frequency-selective  channels.  This  way,  various  reception  schemes  that  are  effective  for 
stationary  frequency-selective  channels  can  be  applied.  In  particular,  we  consider  both  spatial 


multiplexing  schemes  including  RELAX-BLAST  and  Turbo  equalization,  and  space-time  coding 
schemes  including  Alamouti  diversity  techniques. 

Although  Turbo  equalization  significantly  improves  the  detection  performance  compared  to  its 
RELAX-BLAST  counterpart,  its  computational  complexity  is  much  more  expensive  than  RELAX- 
BLAST  mainly  because  to  calculate  the  LMMSE  filter  coefficients,  we  need  to  perform  matrix 
inversion  at  each  time  index.  To  address  this  problem,  we  consider  a  low  complexity  approximation  of 
the  exact  LMMSE  based  Turbo  equalizer,  where  the  approximate  LMMSE  filters  have  constant 
coefficients  over  one  symbol  block.  We  shown  that  Turbo  equalization  employing  such  low 
complexity  equalizer  yields  slightly  degraded  detection  performance  compared  to  that  adapting  a 
directly  implemented  equalizer,  and  this  is  achieved  at  a  computational  complexity  comparable  to 
RELAX-BLAST.  We  hereafter  refer  to  the  Turbo  equalization  scheme  that  employs  a  directly 
implemented  equalizer  as  “Exact-LMMSE-Turbo"  and  its  low  complexity  approximation  as 
“Approximate-LMMSE-Turbo”. 

RESULTS 

1.  WHOI09  Experimental  Results 

The  WHOI09  in-water  experiment  was  conducted  in  December  2009.  The  4  transmit  transducers  with 
source  spacing  up  to  1  m  w'ere  suspended  from  a  vessel  heaving  in  a  14  m  mid-depth  water  column. 
Two  separate  4-hydrophone  arrays  were  deployed  approximately  1  km  and  2  km  away  from  the  source 
array  and  they  will  be  referred  to,  respectively,  as  RBI  and  RB2  receiving  arrays.  Both  RBI  and  RB2 
receiving  arrays  had  0.21  m  spacing  between  adjacent  hydrophones,  and  they  were  mounted  on  an 
anchored  buoy  in  a  mid-water  column  during  the  course  of  data  collection.  The  Doppler  shift  was 
mainly  caused  by  the  relative  motion  between  the  transmitters  and  receivers.  The  carrier  frequency,  the 
sampling  frequency,  and  the  symbol  rate  employed  in  WHO109  experiment  were  30  KHz,  200  KHz, 
and  8  KHz,  respectively. 

GoSLIM  is  employed  to  estimate  the  underlying  CIRs  and  Doppler  frequencies  by  making  use  of  the 
shifted  PeCAN  training  sequences  periodically  allocated  in  the  transmitted  sequence.  Figure  1  shows 
the  CIR  evolution  between  the  active  transmitter  and  the  2nJ  hydrophone  for  all  the  6  epochs  (the 
reason  why  we  consider  the  2nJ  hydrophone  will  be  explained  shortly).  One  observes  that  the  position 
of  the  principal  arrival  shown  in  Figure  1  slowly  shifts  with  time  due  to  the  fact  that  the  transmitter 
array  was  suspended  from  a  vessel  heaving  in  a  water  column  and  the  receiving  hydrophone  arrays 
were  mounted  on  anchored  buoys  during  the  course  of  data  collection.  Figure  2,  on  the  other  hand, 
plots  the  evolution  of  the  Doppler  frequency  for  each  hydrophone.  The  reverberating  Doppler 
frequency  estimates  obtained  from  the  GoSLIM  algorithm  shows  the  presence  of  Doppler  shifts. 

Space-time  cod  ini'  schemes:  Alamouti  coding 

To  investigate  the  performance  of  2Tx-lRx  Alamouti  code,  the  structure  of  the  transmitted  sequences 
is  shown  in  Figure  3.  Two  synchronized  transmitters  were  used.  Each  transmitter  sent  4  data  packets 
and  each  data  packet  was  divided  into  16  blocks  followed  by  a  gap  to  ensure  the  absence  of  inter¬ 
packet  interferences.  In  our  design,  each  payload  segment,  (e.g.,  a  or  b  in  Figure  3),  contains  250 
coded  QPSK  symbols.  Taking  the  segment  a  for  example,  it  is  generated,  as  shown  in  Figure  3,  by 
feeding  250  source  bits  into  a  1/2  rate  convolutional  encoder  followed  by  a  random  interleaver  and 


QPSK  modulation.  Segment  b,  as  well  as  the  segments  in  other  blocks,  is  similarly  generated  but  with 
different  realizations  of  the  source  bits  and  random  interleaver.  The  resulting  uncoded  and  coded  data 
rates  are  7  kbps  and  3.5  kbps,  respectively. 

Since  each  payload  block  contains  a  training  sequence  as  shown  in  Figure  3,  the  reception  scheme  for 
each  block  is  performed  as  follows:  the  receiver  first  conducts  training-directed  channel  estimation 
(note  that  the  GoSLIM  and  GoSLIM-V  algorithms  are  identical  in  this  case  since  only  one  receiver  is 
involved)  by  making  use  of  the  training  sequence,  and  then  detects  the  payload  symbols  contained  in 
the  two  payload  segments  before  and  after  the  training  sequence.  As  a  consequence,  decision-directed 
estimation  is  not  needed.  For  the  sequence  in  Figure  3,  each  block  carries  500  QPSK  payload  symbols. 
Therefore  each  epoch  carries  32  k  symbols  (or  equivalently,  64  k  uncoded  bits  or  32  k  coded  bits)  and 
we  have  3  epochs  for  each  receiving  array.  The  channel  tap  number  is  fixed  at  30  for  all  epochs.  The 
BER  results  at  different  receive  hydrophones,  by  averaging  over  192  k  uncoded  bits  and  96  k  coded 
bits,  are  summarized  in  Table  1  for  each  receiving  array.  One  observes  from  Table  1  that  the  2nd 
hydrophone  of  RBI  array  yields  remarkably  high  BERs.  This  can  be  explained  by  the  fact  that  the 
power  of  the  channel  estimate  at  this  hydrophone  is  persistently  lower  than  that  at  other  hydrophones 
by  almost  an  order  of  magnitude,  which  is  evidenced  in  Figure  I. 

Spatial  multiplexing  schemes^  BLAST_  data 

We  also  designed  4-input  BLAST  data  packet,  which  attained  a  coded  data  rate  of  30  kbps  by 
leveraging  the  MIMO  scheme.  The  payload  sequence  was  divided  into  blocks,  each  of  length  250,  and 
each  block  was  encoded  using  a  'A  convolutional  encoder  with  constraint  length  5.  The  source 
information  is  a  grayscale  Gator  mascot  shown  in  Figure  4(a).  We  made  full  use  of  the  available 
resources  by  transmitting  4  sequences  simultaneously  and  incorporating  the  measurements  acquired 
from  all  the  4  hydrophones,  establishing  a  4x4  MIMO  system.  With  the  RELAX-BLAST  scheme 
employed  in  the  symbol  detection  stage,  we  first  compare  the  resulting  BER  performance  of  the 
GoSLIM  and  GoSLIM-V  algorithms  in  the  channel  estimation  stage.  When  GoSLIM  is  employed, 
Figures  4(b),  4(c),  and  4(d)  show  the  recovered  grayscale  mascots  for  epochs  “195600”,  “195730”,  and 
“195860”,  respectively.  The  corresponding  coded  BERs  and  the  computational  time  consumed  at  the 
channel  estimation  stage  on  an  ordinary  workstation  (Intel  Xeon  E5506  processor  2.13G  Hz,  12GB 
RAM,  Windows  7  64-bit,  and  MATLAB  R2010b)  are  listed  in  the  first  row  of  Tables  2  and  3, 
respectively.  In  comparison,  when  GoSLIM-V  is  employed,  the  recovered  grayscale  mascots  are 
shown  in  Figures  4(e),  4(f),  and  4(g),  with  the  corresponding  BERs  and  computational  complexities 
listed  in  the  second  row  of  Tables  2  and  3,  respectively.  One  observes  from  Table  2  that  GoSLIM-V 
yields  slightly  better  BER  results  than  its  GoSLIM  counterpart.  Moreover,  Table  3  demonstrates  that 
GoSLIM-V  is  about  4  times  faster  than  GoSLIM.  Due  to  the  advantage  of  GoSLIM-V  over  the  original 
GoSLIM  algorithm,  we  hereafter  employ  GoSLIM-V  in  the  channel  estimation  stage  and  compare  the 
RELAX-BLAST  scheme  with  the  Turbo  equalization  counterpart  in  the  symbol  detection  stage. 

We  remark  that  the  three  grayscale  mascots  carried  by  epochs  “195600”,  “195730”,  and  “195860”, 
respectively,  can  be  perfectly  recovered  by  employing  either  type  of  the  Turbo  equalization  scheme. 
Moreover,  Table  4  lists  the  ratio  of  the  time  consumed  to  process  a  packet  to  the  duration  of  the  same 
packet  during  transmission,  obtained  using  Exact-LMMSE-Turbo,  Approximate-LMMSE-Turbo,  and 
the  RELAX-BLAST  scheme,  respectively.  The  results  are  obtained  by  applying  3  iterations  for  all  of 
the  three  types  of  detection  schemes  considered. 


2.  ACOMM10  Experimental  Results 


The  ACOMMIO  experiment  took  place  in  July  2010  in  the  Mid-Atlantic  Bight  on  the  continental  shelf 
off  the  coast  of  New  Jersey  in  an  area  with  water  depth  of  78  m.  The  transmitter  array  consisted  of  12 
transducers  with  0.8  in  spacing  between  adjacent  elements.  In  our  design,  only  4  transducers  were 
activated.  The  receiving  array  was  composed  of  8  hydrophones.  The  spacing  between  the  adjacent 
hydrophones  was  2.06  m  except  the  first  element,  which  was  spaced  4  m  above  the  second  element. 
Both  the  transmit  array  and  the  receiving  array  were  mounted  on  anchored  buoys  and  were  deployed 
approximately  3  km  away  from  each  other.  The  Doppler  shift  was  mainly  caused  by  the  relative 
motion  between  the  transmitters  and  receivers.  The  carrier  frequency,  the  sampling  frequency,  and  the 
symbol  rate  employed  in  the  ACOMMIO  experiment  were  20  kHz,  80  kHz,  and  4  kHz,  respectively. 

The  BLAST  data  package  of  20  s  in  duration  was  transmitted  multiple  times  in  the  ACOMMIO 
experiment  and  was  recorded  by  the  receiving  array.  A  total  of  89  epochs  were  available  and  they  are 
referred  to  as  "M1MO01”— “M1M089”,  respectively.  Figures  5  and  6  show,  respectively,  the 
evolutions  of  CIR  and  Doppler  frequencies  estimated  by  GoSLIM.  The  position  shift  of  the  principal 
arrival  with  respect  to  time  shown  in  Figure  5  and  the  reverberating  Doppler  frequency  estimates 
shown  in  Figure  6  suggest  the  presence  of  the  Doppler  shifts. 

Sjnitujl  multiplexine  schemes:  BLASTS  data 

Figure  7  shows  the  source  information  contained  in  the  transmitted  package.  Each  package  consists  of 
7  packets.  The  first  4  packets  convey  4  grayscale  Gator  mascots  and  the  last  3  packets  combined  form 
a  colored  mascot.  The  RGB  components  of  the  colored  image  are  transmitted  in  the  5lh ,  6,h,  and  7,h 
packets,  respectively.  Guard  intervals  between  adjacent  packets  are  used  to  prevent  inter-packet 
interference.  Each  pixel  of  the  grayscale  image  is  represented  by  5  bits,  corresponding  to  32  different 
intensities  (e.g.,  pure  white  and  pure  dark  pixels  are  represented  by  11111  and  00000,  respectively). 
The  64-pixel  by  1 00-pixel  grayscale  mascot  image,  as  a  consequence,  is  represented  by  a  total  of  32  k 
source  bits.  (Accordingly,  a  colored  mascot  image  is  represented  by  96  k  bits.)  The  contrast  of  the 
grayscale  image,  as  well  as  the  hue  of  the  colored  image,  has  been  carefully  adjusted  so  that  the  image 
carries  approximately  equal  numbers  of  1  ’s  and  0's.  The  source  information  obtained  from  the  Gator 
mascots  was  divided  into  blocks,  each  of  length  250,  and  each  block  was  encode  using  a  Vi 
convolutional  encoder  with  constraint  length  5.  By  transmitting  N  =  4  sequences  simultaneously  and 
incorporating  the  measurements  acquired  from  all  of  the  A/  =  8  receiver  elements  for  analysis,  we 
established  a  4x8  MIMO  UAC  system,  which  leads  to  a  net  coded  data  rate  of  15  kbps. 

We  first  assess  the  impact  of  the  channel  estimation  algorithm  (GoSLIM  or  GoSLIM-V)  on  the 
resulting  detection  performance.  With  the  estimated  CIRs  and  Doppler  frequency  (frequencies),  after 
phase  compensation,  the  payload  symbols  are  detected  using  the  RELAX-BLAST  scheme.  There  is  a 
total  of  623  packets  (recall  that  we  have  89  epochs  and  each  consists  of  7  packets)  available  and  we 
deem  a  packet  to  be  successfully  detected  if  its  coded  BER  is  less  than  0.1.  When  GoSLIM  is 
employed  as  the  channel  estimation  algorithm,  we  have  succeeded  in  tracking  the  entire  32  payload 
blocks  for  594  packets.  A  coded  BER  of  5.1x1 0-1  is  achieved  after  averaging  over  the  594  successful 
packets.  In  comparison,  when  GoSLIM-V  is  used,  599  packets  are  successfully  retrieved  with  an 
average  coded  BER  of  4.9  xlO-3.  We  use  the  1st  packet  (i.e.,  the  one  conveys  a  grayscale  gator  mascot; 
see  Figure  7)  of  three  epochs,  namely  “MIMO02”,  “M1MO06",  and  “M1M018”,  to  compare  the 
impact  of  the  two  channel  estimation  algorithms  on  the  resulting  BER  performance  and  computation 
complexity.  When  GoSLIM  is  employed,  Figures  8(a),  8(b),  and  8(c)  show  the  recovered  grayscale 


mascots  for  “MIMO02”,  “MIMO06”,  and  “MIM018”,  respectively.  The  corresponding  coded  BERs 
and  the  computational  time  consumed  at  the  channel  estimation  stage  on  an  ordinary  workstation  (Intel 
Xeon  E5506  processor  2.13G  Hz,  12GB  RAM,  Windows  7  64-bit,  and  MATLAB  R2010b)  are  listed 
in  the  first  row  of  Tables  5  and  6,  respectively.  In  comparison,  when  GoSLIM-V  is  employed,  the 
recovered  grayscale  mascots  are  shown  in  Figures  8(d),  8(e),  and  8(f),  with  the  corresponding  BERs 
and  computational  complexities  listed  in  the  second  row  of  Tables  5  and  6,  respectively.  One  observes 
from  Table  5  that  GoSLIM-V  yields  slightly  better  BER  results  than  its  GoSLIM  counterpart. 
Moreover,  Table  6  demonstrates  that  GoSLIM-V  is  about  4  times  faster  than  GoSLIM.  Due  to  the 
advantage  of  GoSLIM-V  over  the  original  GoSLIM  algorithm,  we  hereafter  employ  GoSLIM-V  in  the 
channel  estimation  stage  and  compare  the  RELAX-BLAST  scheme  with  the  Turbo  equalization 
counterpart  in  the  symbol  detection  stage. 

Tabic  7  lists  the  successfully  detected  packet  percentage,  the  error-free  packet  percentage,  the  coded 
BER  averaged  over  the  successful  packets,  and  the  time  ratio  of  the  time  consumed  to  process  a  packet 
to  the  duration  of  the  same  packet  during  transmission,  obtained  using  Exact-LMMSE-Turbo, 
Approximate-LMMSE-Turbo,  and  the  RELAX-BLAST  scheme,  respectively.  The  results  are  obtained 
by  apply  ing  3  iterations  for  all  of  the  three  types  of  detection  schemes  considered.  One  observes  from 
Table  7  that  1)  BER-wise,  both  Exact-LMMSE-Turbo  and  Approximate-LMMSE-Turbo  outperform 
RELAX-BLAST  significantly,  2)  Approximate-LMMSE-Turbo  greatly  reduces  the  complexity  at  the 
cost  of  slight  BER  performance  degradation  compared  to  Exact-LMMSE-Turbo,  and  3)  compared  to 
RELAX-BLAST,  Approximate-LMMSE-Turbo  improves  the  BER  performance  by  two  orders  of 
magnitude  without  significantly  increasing  the  computational  complexities. 


3.  MACE10  Experimental  Results 

The  MACEIO  in-water  experiment  was  conducted  off  the  coast  of  Martha's  Vineyard,  MA  in  June 
2010.  A  source  array  consisting  of  4  transducers  was  vertically  deployed  at  a  depth  of  80  m,  and  towed 
by  a  vessel.  At  the  receiver  side,  a  12-element  hydrophone  array  was  mounted  on  a  buoy.  The  vessel 
moved  from  the  minimum  range  of  500  m  away  from  the  receiving  array  outbound  to  the  maximum 
range  of  4000  m  and  then  inbound  back  to  the  minimum  range.  The  carrier  frequency,  the  sampling 
frequency,  and  the  symbol  rate  employed  in  the  MACEIO  experiment  were  13  kHz,  39.0625  kHz,  and 
3.90625  kHz,  respectively.  MACEIO  is  a  mobile  UAC  experiment,  where  Doppler  induced  temporal 
scaling  effects  (stretching  or  compression)  are  present.  Temporal  resampling  is  thus  implemented  to 
convert  the  Doppler  scaling  effects  to  Doppler  frequency  shifts.  The  channel  model  remained  after  the 
temporal  resampling  procedure  is  similar  to  that  in  the  WHIO09  and  ACOMM10  experiments.  We 
compare  the  GoSLIM  and  GoSLIM-V  algorithms  in  the  channel  estimation  stage  and  compare  the 
RELAX-BLAST  scheme  and  the  LMMSE  Turbo  equalization  technique  (including  its  low  complexity 
approximation)  in  the  symbol  detection  stage. 

The  structure  of  a  transmitted  data  package  is  shown  in  Figure  9,  which  is  similar  to  that  of  the 
ACOMM10  experiment.  Each  package  consists  of  4  packets  which  are  designed  to  test  the  MIMO 
BLAST  scheme.  The  1st  packet  conveys  a  grayscale  Gator  mascot  and  the  subsequent  3  packets 
combined  form  a  colored  mascot.  The  RGB  components  of  the  colored  image  were  transmitted  in  the 
2nd,  3rd,  and  4th  packets,  respectively.  Each  packet  is  constructed  as  follows  (see  Figure  9):  time- 


marking  sequences  are  placed  at  the  beginning  of  each  packet  to  facilitate  the  temporal  resampling 
procedure;  two  guard  intervals  (each  containing  500  silent  symbols)  are  placed,  respectively,  before 
and  after  the  segments  containing  the  payload  symbols  and  training  sequences.  The  payload  symbols 
contain  the  information  of  the  Gator  mascot  image.  By  transmitting  N  -4  sequences  simultaneously 
and  incorporating  the  measurements  acquired  from  all  of  the  M  =  12  receiver  elements  for  analysis, 
we  established  a  4x12  M1MO  UAC  system,  which  leads  to  a  net  coded  data  rate  of  1 1.7  kbps.  The 
data  package  was  transmitted  periodically  and  recorded  by  the  receiver  array.  A  total  of  120  epochs 
were  available  and  they  are  referred  to  as  “E001”— “E120”,  respectively. 

The  packet  is  resampled  using  the  Doppler  scaling  factor  estimated  by  taking  advantage  of  the  time¬ 
marking  sequences.  To  assess  the  performance  of  the  resampling  process,  the  CIR  and  Doppler 
frequency  evolutions  obtained  by  GoSLIM  before  we  resample  the  2nJ  packet  of  epoch  ‘‘E002”  are 
shown  in  Figures  10(a)  and  10(c),  respectively.  In  comparison,  Figures  10(b)  and  10(d)  demonstrate 
the  corresponding  CIR  and  Doppler  frequency  evolutions  obtained  after  resampling  the  packet, 
respectively.  We  can  see  from  Figure  10  that  the  temporal  resampling  process  successfully  reduces  the 
Doppler  scaling  effects  to  Doppler  frequency  shifts.  Moreover,  it  is  interesting  to  look  at  Figure  1 1 
where  the  vessel  speed  estimated  during  the  resampling  stage  is  plotted  on  top  of  the  GPS  reference 
information  provided  by  WHOl  (the  GPS  device  was  equipped  on  the  moving  vessel).  The  good 
agreement  between  these  two  curves  verifies  the  effectiveness  of  the  resampling  procedure  we  employ. 
The  analysis  presented  hereafter  is  based  on  the  resampled  measurements. 

Spatial  multiplexing,  schemes^  RLAST_  data 

We  choose  the  lsl  packet  of  epoch  “E018”  to  verify  the  channel  model  behind  GoSLIM-V  (other 
epochs  give  similar  observations).  The  evolution  of  the  Doppler  frequencies  produced  by  GoSLIM  for 
all  the  12  receive  hydrophones  are  plotted  superimposed  in  Figure  12  along  with  the  evolution  of  the 
Doppler  frequency  obtained  by  GoSLIM-V.  One  observes  that  the  curves  show  good  agreement  with 
each  other,  which  verities  the  validity  of  the  key  assumption  for  GoSLIM-V  that  different  receivers 
experience  the  same  Doppler  frequency.  We  proceed  to  assess  the  impact  of  the  channel  estimation 
algorithm  (GoSLIM  or  GoSLIM-V)  on  the  resulting  detection  performance.  With  the  estimated  CIRs 
and  Doppler  frequency  (frequencies),  after  phase  compensation,  the  payload  symbols  are  detected 
using  the  RELAX-BLAST  scheme. 

There  is  a  total  of  480  packets  available  and  we  deem  a  packet  to  be  successfully  detected  if  its  coded 
BER  is  less  than  0.1.  When  GoSLIM  is  employed  as  the  channel  estimation  algorithm,  we  have 
succeeded  in  tracking  the  entire  32  payload  blocks  for  391  packets.  A  coded  BER  of  1.7x1  O'2  is 
achieved  after  averaging  over  the  391  successful  packets.  In  comparison,  when  GoSLIM-V  is  used, 
396  packets  are  successfully  retrieved  with  an  average  coded  BER  of  1.6x1  O'2.  We  use  the  lsl  packet 
(i.e.,  the  one  conveys  a  grayscale  gator  mascot;  see  Figure  9)  of  three  epochs,  namely  “E013”,  “E016”, 
and  “E018”,  to  assess  the  impact  of  the  two  channel  estimation  algorithms  on  the  resulting  BER 
performance  and  computation  complexity.  When  GoSLIM  is  employed,  Figures  13(a),  13(b),  and  13(c) 
show  the  recovered  grayscale  mascots  for  '‘E013’\  “E016”,  and  “E018’\  respectively.  The 
corresponding  coded  BERs  and  the  computational  time  consumed  at  the  channel  estimation  stage  on  an 
ordinary  workstation  (Intel  Xeon  E5506  processor  2.13G  Hz,  12GB  RAM,  Windows  7  64-bit,  and 
MATLAB  R2010b)  are  listed  in  the  first  row  of  Tables  8  and  9,  respectively.  In  comparison,  when 
GoSLIM-V  is  employed,  the  recovered  grayscale  mascots  are  shown  in  Figures  13(d),  13(e),  and  13(f), 
with  the  corresponding  BERs  and  computational  complexities  listed  in  the  second  row  of  Tables  8  and 


9,  respectively.  One  observes  from  Table  8  that  GoSLIM-V  yields  slightly  better  BER  results  than  its 
GoSLIM  counterpart.  Moreover,  Table  9  demonstrates  that  GoSLIM-V  is  about  4  times  faster  than 
GoSLIM.  Due  to  the  advantage  of  GoSLIM-V  over  the  original  GoSLIM  algorithm,  we  hereafter 
employ  GoSLIM-V  in  the  channel  estimation  stage  and  compare  the  RELAX-BLAST  scheme  with  the 
Turbo  equalization  counterpart  in  the  symbol  detection  stage. 

Table  10  lists  the  successfully  detected  packet  percentage,  the  error-free  packet  percentage,  the  coded 
BER  averaged  over  the  successful  packets,  and  the  time  ratio  of  the  time  consumed  to  process  a  packet 
to  the  duration  of  the  same  packet  during  transmission,  obtained  using  Exact-LMMSE-Turbo, 
Approximate-LMMSE-Turbo,  and  the  RELAX-BLAST  scheme,  respectively.  The  results  are  obtained 
by  applying  3  iterations  for  all  of  the  three  types  of  detection  schemes  considered.  Similar  observations 
can  be  obtained  from  fable  10  as  from  Table  7  in  the  ACOMMIO  experimental  results.  Moreover,  we 
analyze  epoch  “E054”  that  leads  to  perfect  recovery  of  both  the  grayscale  and  colored  mascots  (sec 
Figure  14(b)  and  14(d))  using  either  Exact-LMMSE-Turbo  or  Approximate-LMMSE-Turbo.  In 
comparison,  the  grayscale  and  colored  mascots  recovered  from  epoch  “E054’’  using  RELAX-BLAST 
are  shown  in  Figures  14(a)  and  14(c),  respectively,  with  the  corresponding  coded  BERs  being 
1.8x10-'  and  1.6x10-'. 


IMPACT/APPLICATIONS 

The  natural  bandwidth  limitations  of  coherent  underwater  acoustic  channel  suggest  a  technical 
breakthrough.  MIMO  signal  processing  is  a  promising  bandwidth  efficient  method  to  high  data  rate 
and  high  quality  serv  ices.  Our  promising  results  are  expected  to  favorably  impact  high-rate  long-range 
MIMO-UAC  designs. 
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Figure  1:  CIK  evolution  over  approximately  a  9  second  period  at  the  2nd  receiver  for  all  the  6  epochs . 
[graph:  One  observes  from  the  GoSL/M  CIR  estimates  that  the  position  of  the  principal  arrival 
shifts  with  time .  This  observation  is  in  line  with  the  fact  that  the  transmitter  array  was  suspended 
from  a  vessel  heaving  in  a  14  m  mid-depth  water  column  and  the  receiving  hydrophone  arrays  were 
mounted  on  anchored  buoys  during  the  course  of  data  collection  in  WHO109  in-water  experiment] 
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Figure  2:  Evolution  of  the  estimated  Doppler  frequencies  at  each  receiver  for  all  the  6  epochs, 
[graph:  The  reverberating  Doppler  frequency  estimates  obtained  from  the  GoSLIM  algorithm 
provide  evidence  on  the  presence  of  Doppler  shifts.] 
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Figure  3:  The  structure  of  the  transmitted  sequence  using  2Tx~lRx  Aiamouti  scheme . 
[graph:  Each  transmitter  sends  4  data  packets  and  each  packet  is  divided  into  16  blocks.  The 
payload  symbols  contained  in  each  block  are  constructed  using  Aiamouti  diversity  scheme.] 


RBI 

RB2 

un coded  BER 

coiieil  BER 

uncoded  BER 

coded  BER 

Rxl 

2.7083  x  l()-“ 

0 

1.8750  x  10 

0 

Rx2 

0.0180 

0.0021 

1.2500  x  lO"4 

0 

Rx3 

1.6146  x  10~* 

0 

6.2498  x  10“5 

0 

Rx4 

3.2291  x  IQ”4 

0 

1.1458  x  10~* 

0 

Table  1:  Coded  and  uncoded  BER  performance  of  GoSLIM  coupled  with  Aiamouti  diversity  scheme 

for  2Tx-lRx  system. 

[graph:  GoSLIM  coupled  with  Aiamouti  diversity  scheme  yields  satisfactory  detection  performance 
except  for  the  2nd  hydrophone  oj  RBI  receiving  array.  The  reason  is  that  the  estimated  channel 
power  of  this  hydrophone  is  persistently  lower  than  those  of  other  hydrophones  by  almost  an  order 

of  magnitude ,  see  Figure  L] 
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Figure  4:  (a)  The  grayscale  Gator  mascot  transmitted  (truth),  (b)  Mascot  recovered  from  Epoch 
“195600-RB2”,  (c)  Mascot  recovered  from  Epoch  “ 195730-RB2 ”.  (d)  Mascot  recovered  from  Epoch 
“195H60-RB2”.  (e)  Mascot  recovered  from  Epoch  “I9S600-RB2  (f)  Mascot  recovered  from  Epoch 
*‘195730-RB2”.  (g)  Mascot  recovered  from  Epoch  “195860-RB2”.  (b)-(d)  are  obtained  by  GoSLIM. 

(e)-(g)  are  obtained  by  GoSLIM-V. 


1 95600 

195730 

195860 

GoSLIM 

2.0  x  in-3 

0.7  x  10-  1 

7.2  x  10— 1 

GoSLIM-V 

1.8  x  10-3 

8.1  x  10— 1 

5.0  x  10— 1 

Table  2:  Coded  BER  results  obtained  by  GoSLIM  and  GoSLIM-  V,  respectively. 


Time  {sec)  Time  (sec) 


195600 

195730 

195860 

GoSLIM  (15  iterations) 

38.6 

39.6 

40.1 

GoSLIM-V  (15  iterations) 

10.5 

10.7 

10.8 

Table  3:  Complexity  comparison  (in  second)  between  GoSLIM  and  GoSLIM-V. 


195600 

195730 

195S60 

Exact-LMMSE-Turbo  (3  iterations) 

512.0 

513.3 

514.0 

Approximate-LMMSE-Turbo  (3  iterations) 

21.9 

22.3 

22.5 

1  '  . 

RELAX-BLAST  (3  iterations) 

19.8 

20.0 

20.5 

Table  4:  The  time  ratio  summary  of  Exact-LMMSE-Turbo,  Approximate-  LMMSE-Turbo,  and 

RELAX-BLAST. 
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Figure  5:  CIR  evolutions  obtained  by  GoSLIM  between  the  four  active  transmitters  and  one 

hydrophone  for  Epoch  “MIM028” . 
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Figure  6:  Evolutions  of  the  Doppler  frequencies  estimated  by  GoSLIM  at  each  receiver  for  Epoch 

“ MIMQ28 
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Figure  7:  Each  package  transmitted  in  the  ACOMM 10  experiment  contains  4  grayscale  Gator 
mascot  images  and  1  colored  image.  The  colored  image  is  decomposed  into  RGB  components. 
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Figure  8:  (a)  Grayscale  mascot  recovered  from  epoch  “MIMO02”.  (b)  Grayscale  mascot  recovered 
from  epoch  “MIMO06”.  (c)  Grayscale  mascot  recovered  from  epoch  “MIM018”.  (d)  Grayscale 
mascot  recovered  from  epoch  “MIMO02”.  (e)  Grayscale  mascot  recovered  from  epoch  “MIMO06”. 
(f)  Grayscale  mascot  recovered  from  epoch  “MIM018”.  (a)-(c)  are  obtained  by  GoSLIM.  (d)-(f)  are 

obtained  by  GoSLIM-V. 


MIMO02 

MIMO06 

MIM01S 

GoSLIM 

(>.2  x  Hr3 

7.8  x  10“  1 

:L9  x  10“ 3 

GoSLIM-V 

5.4  x  1<)-3 

G.3  x  10—' 

1.2  x  10“ 3 

Table  5:  Coded  BER  results  obtained  by  GoSLIM  and  GoSLIM-V,  respectively. 


MIMO02 

MIMO06 

MIM018 

GoSLIM  (15  iterations) 

82.0 

81.3 

83.1 

GoSLIM-V  (15  iterations) 

21.4 

20.5 

21.2 

Table  6:  Complexity  comparison  (in  second)  between  GoSLIM  and  GoSLIM-V. 


successful  packet  percentage  (*•) 

crror-irec  packet  percentage  (*•) 

average  coded  BER 

tune  ratio 

Exact -LMMSE-Turbo  (5  itcmoiis) 

100 

100 

0 

442.7 

Approximate- LMMSE-Turbo  (3  derations) 

100 

9  T9 

1.6  x  10  ”5 

152 

RELAX-BLAST  (3  iterations  1 

95  3 

S.7 

5.1  >  10  3 

14.5 

Table  7:  The  performance  summary  of  Exact-LMMSE-Turbo,  Approximate-  LMMSE-Turbo,  and 

RELAX-BLAST. 
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Figure  9:  The  structure  of  the  package  used  in  the  MACE  10  experiment. 
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Figure  10:  (a)  CIR  evolution  of  Epoch  “E002”  before  resampling,  (b)  CIR  evolution  of  Epoch 
“E002”  after  resampling.  (c)  Doppler  frequency  evolution  of  Epoch  “E002”  before  resampling,  (d) 
Doppler  frequency  evolution  of  Epoch  “E002”  after  resampling. 


Figure  11:  The  relative  speed  between  the  transmitter  and  receiver  array  given  by  GPS  and 
estimated  during  the  temporal  resampling  stage  (courtesy  of  Milica  Stojauovic's  group). 


Figure  12:  Doppler  frequency  evolution  of  the  first  packet  in  epoch  “ED  18"  obtained  by  GoSLIM 

and  GoSLIM-V,  respectively. 


E013 

E016 

E018 

GoSLIM 

3.6  x  10~2 

3.6  x  10-3 

1.6  x  10~4 

GoSLIM-V 

2.3  x  10-2 

3.1  x  10"3 

9.4  x  IQ"5 

Table  8:  Coded  RER  results  obtained  by  GoSLIM  and  GoSLIM-V,  respectively. 


E013 

E016 

E018 

GoSLIM  (15  iterations) 

122.8 

121  1 

123.0 

GoSLIM-V  (15  iterations) 

32.1 

31.8 

32.3 

Table  9:  Complexity  comparison  (in  second)  between  GoSLIM  and  GoSLIM-V. 
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Figure  13:  (a)  Grayscale  mascot  recovered  from  epoch  “ EG  13 (h)  Grayscale  mascot  recovered 
from  epoch  “11016".  (c)  Grayscale  mascot  recovered  from  epoch  “E018’\  (d)  Grayscale  mascot 
recovered  from  epoch  “E013  (e)  Grayscale  mascot  recovered  from  epoch  “E016”.  (f)  Grayscale 

mascot  recovered  from  epoch  “E018".  (a)-(c)  are  obtained  by  GoSLIM.  (d)-(f)  are  obtained  by 

GoSLIM-V. 
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l(K) 
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S.l  x  10~4 

17.2 
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4  S 

1.6  x  I(T2 
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Table  10:  The  performance  summary  of  Exact-LM MS  E-Turbo,  Approximate-  LMMSE-Turbo,  and 

RELAX-BLAST. 


(c)  (d) 

Figure  14:  (a)  Grayscale  mascot  recovered  from  epoch  “E054”  using  RELAX-BLAST.  (b) 
Grayscale  mascot  recovered  from  epoch  “E054”  using  Turbo  equalization,  (c)  Colored  mascot 
recovered  from  epoch  “ E054 ”  using  RELAX-BLAST.  (d)  Colored  mascot  recovered  from  epoch 

“E054”  using  Turbo  equalization. 
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